s17adf
s17adf
© Numerical Algorithms Group, 2002.
Purpose
S17ADF Bessel function Y (x)
1
Synopsis
[fun_value,ifail] = s17adf(x<,ifail>)
Description
This routine evaluates an approximation to the Bessel Function of
the second kind Y (x).
1
Note: Y (x) is undefined for x<=0 and the routine will fail for
1
such arguments.
2
For x near zero, Y (x)~=- -----. This approximation is used when
1 (pi)x
x is sufficiently small for the result to be correct to machine
precision. For extremely small x, there is a danger of overflow
2
in calculating - ----- and for such arguments the routine will
(pi)x
fail.
For very large x, it becomes impossible to provide results with
any reasonable accuracy, hence the routine fails. Such arguments
contain insufficient information to determine the phase of
______
/ 2
oscillation of Y (x), only the amplitude, / -----, can be
1 \/ (pi)x
determined and this is returned on soft failure. The range for
which this occurs is roughly related to machine precision; the
routine will fail if x>~1/machine precision.
Parameters
s17adf
Required Input Arguments:
x real
Optional Input Arguments: <Default>
ifail integer -1
Output Arguments:
ifail integer